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^ ■ We investigate the statistical properties of single phase slip events observed 

, in vortex nucleation experiments in ^He by the groups in Berkeley and Paris. 

ly-) [ From the cumulative distribution function of the events we calculate the slip 

■ rate as a function of flow velocity. The critical velocity is defined as the mean 

velocity and its statistical width as the standard deviation. From the slip rate 
and from the observed linear temperature dependence of the critical velocity 
we obtain the energy barrier for vortex nucleation which is a quadratic func- 
tion of the flow velocity. A comparison with the statistical properties of the 
laminar to turbulent transition in the flow around an oscillating sphere shows 
'"O ' strikingly different behaviour. 

PACS numbers: 67.40. Vs, 67.40. Hf. 

^ . 1. Introduction 

Nucleation of individual quantized vortices in the flow of superfluid 4 He 
through micro-apertures has been observed since 1985J1H3'@ In these beauti- 
ful and meanwhile famous experiments the 2tt phase slip of the superfluid 
wavefunction is detected which occurs when a vortex crosses the stream lines 
of the aperture. The nucleation event was shown to be a stochastic process 
depending on the flow velocity and the temperature. Consequently, a sta- 
tistical analysis of the events is performed. Quantities like the cumulative 
distribution function, slip rate, critical velocity and its statistical width are 
introduced and analyzed in order to gain insight into the physics of the nu- 
cleation process. One of the most interesting properties is the energy barrier 
which a vortex has to overcome during nucleation either by thermal activa- 
tion or by quantum tunneling. The dependence of the barrier on the flow 
velocity is of particular importance for developing a physical model for the 
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creation of vortices. 

In our present work we are investigating the statistical properties of the 
nucleation events in a different way. The motivation comes from our own 
experiments on the transition from laminar flow to turbulence around an 
oscillating micro-sphere in superfluid 4 He at millikelvin temperaturesE The 
onset of turbulence in that case is a stochastic process which can be de- 
scribed appropriately by means of reliability theory. With this tool at hand 
we are re-analyzing the statistical properties of the single vortex nucleation 
events. Our method is simple, direct and transparent. It avoids mathemati- 
cal shortcomings of the previous analysis by the authors. From the rigorous 
relation between cumulative distribution function and slip rate we find the 
latter as a function of flow velocity v. From this result we then infer the 
energy barrier E{v) which is consistent with the observed linear temperature 
dependence of the critical flow velocity v c . We are analyzing both ac and dc 
flow experiments and find similarities as well as differences. Finally, we com- 
pare our results of the phase slip experiments with the statistical analysis of 
our own experiment. 



2. Basics of Reliability Theory 

In this Section we present the basic tools of reliability theory! needed for 
the statistical analysis of the phase slip events. The cumulative distribution 
function (CDF) of the probability of observing a slip below a given velocity 
amplitude v is P(v), which is obtained from recording a time series of slip 
events. In order to employ the terminology of reliability theory we consider 
"slip" as "failure" and "no-slip" as "reliability". The reliability function 

R(v) = 1 - P(v) 

is determined by the failure rate A(v): 



R(v) = exp J -J A{v')dv' 



from which follows 



, . d InR 

Hv) = -r- ■ (1) 

dv 

We note that A has the dimension of 1 /velocity, i.e., s/m. It can be simply 
transformed to units of 1/s by the relation 

. I * , . . dv 
A(t)=A(,)-, 
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where dv/dt is the change of the flow velocity amplitude with time. Be- 
cause in the experiments the amplitude increases linearly, dv/dt is a con- 
stant during the time series and depends only on drive and dissipation. The 
probability density function (PDF) f(y) is given by 

, dP dR 

f{v) = *7 = " *7 • 

From the experiment f{v) is obtained if histograms are evaluated. The 
failure rate may be calculated directly from f{v) by using 

AW- '<"> 



//(»')*' 

v 

The critical velocity v c of the flow may be defined as the mean velocity 



'Jc = {v) = J v ■ f(v)dv . 



Alternatively, the median velocity v m may be used which is given by P(v m ) = 
R(v m ) = 1/2 or by 

fin 

f(v)dv = - . 
o 

In the following it is irrelevant which choice is used. Finally, because the 
transition of P from to 1 occurs over a finite velocity interval, a "width" 
Av c may be defined. We use the standard deviation 



v c = \l (v 2 ) - {vY 



A 

i.e., the square-root of the variance. This choice is numerically much safer 
than the authors' definition of the inverse slope of P(v) at P = 1/2 which 
is given by l//(f m ) and therefore depends only on one point on the PDF 
which may easily be subject to errors. 



3. Data Analysis 

In order to apply the above theory to a recorded time series we can start 
either from the CDF or from the PDF. In the published literature these data 
can be found in two places, namely in the recent review by the Paris groupcl 
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Fig. 1. Reliability function R(v) from the Berkeley PDF at 0.46 K (bottom 
axis) and from the Berkeley CDF at 0.79 K (top axis) (see Ref.|3|). Note 
the temperature independent shape. The overlap is obtained by shifting the 
velocity scales by 0.84 m/s with respect to each other. 
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Fig. 2. Slip rate A(v) obtained from the data in Fig.g: a) 0.46 K, b) 0.79 K. 
Note the exponential variation having the same slope for both temperatures. 



where data at 12 mK are shown, and in a publication by the Berkeley grou 
where a PDF at 0.46 K and a CDF at 0.79 K are displayed. We have scanned 
these data and digitized them for analysis. 

A first very interesting observation is the temperature independence of 
the shape of the Berkeley CDFs. Both curves overlap if the velocity scale is 
shifted, see Fig. [j]. This shift is due to a temperature dependence of v c (or 
v m ) which is linearB At 0.46 K and at 0.79 K we find v c = 5.07 m/s and 
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4.23 m/s, respectively. The median velocity v m is larger by only 0.01 m/s 
at either temperature. The constant shape implies a constant width Av c 
= 0.0611 m/s. The slip rate A(v) is obtained by numerical differentiation, 
see Eq. (||). The results for both temperatures are shown in Fig. ||. An 
exponential increase of A is obvious with a temperature independent slope b : 

A(v,T) =a(T) -exp(b-v) , 

where b = 20.5 s/m (1/6 = 0.0488 m/s) and a(T) is a temperature depen- 
dent coefficient. Inspection of a(T) gives the following surprisingly simple 
relation: 

a(T) = T v • exp(-6 • v c (T)), 
where T v = 14.3 s/m. Therefore, we have finally 

A(v,T) = rVexp{&>- Uc (T))}, (2) 

and at all v = v c (T) we obtain A(i> c ) = T v . The critical velocity is the 
"working point": because of the rapid variation of A(v) data can be taken 
only in the vicinity of v c . The time resolution of the experiment is limited 
by the half period of the oscillating flow because the amplitude of each 
half period is recorded. Therefore, the highest measurable slip rate in the 
Berkeley ac experiments is given by A(t) =2/, where / = 40.7 Hz is the 
oscillation frequency, or A(v) = A(t)/ (dv/dt) = 581 s/m where we used the 
value dv/dt = 0.14 m/s 2 obtained from Fig.||] of Ref. ||. In fact, the original 
data appear to extend up to about this value. Because of the digitizing 
procedure our data are limited to about 200 s/m. 

Performing the same analysis with the Paris CDf! we obtain the same 
exponential behaviour of A(v). Because in that work the velocity is given 
in dimensionless units, we cannot compare the absolute values directly with 
the above results from the Berkeley data. We find: v c = 59.64, Av c = 0.195, 
b = 6.23, T v = 3.4. The highest slip rate (measured at v = 60) is A(60) = 15 
and from dv/dt ~ 2.1 s _1 (from Fig.|| of Ref.^) we get the maximum rate 
Ait) = 32 s _1 which corresponds to one slip per half period (which is 31.8 
ms).i 

From the expression of the slip rate (see Eq.@) we easily obtain an 
analytical form of R(v) and f(v). Neglecting the small difference between 
v c and v m we find from R(v m ) = 1/2 the following results: f(v m ) = T v /2 
and T v = b ■ In 2. Both the Berkeley and the Paris data are consistent with 
this. It is quite plausible that only two parameters are left to determine 
the properties of the CDF or the PDF: the critical velocity v c (T) , which 
determines the position on the velocity axis, and the width Ai> c , which can 
be expressed either by T v or by b. 
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We now draw conclusions from the temperature independence of the 
slope b. Because the ac data can be taken only in the vicinity of v c it is 
reasonable to compare Eq. (g) with a Taylor expansion of the energy barrier 
E(v) for vortex nucleation which has to be overcome by thermal activation:! 

E(v)-E(v c ) = (^-) -(v-v c ). 



dv /, 

Comparing —(E(v) — E(v c ))/T with the exponent in Eq. ([|) we find 



1 fdE 

dv 



vc(T) 



which according to the experiment must be independent of T. This implies 
a differential equation 

'dE\ 

-) =-bT 

Vc(T) 



dv J. 



where T is the inverse function of v c (T) = v c q(1 — T/Tq) which is obtained 
from the experiment. Integration then gives 

E(v)=E -T bv + ^v 2 

for all v = v c (T), i.e., for < v < v c q. Choosing E(v c o) = fixes the 
integration constant to the value Eq = bv c Q Tq/2 and gives 

E(v) = E (l-—) 2 . (3) 

Inserting 6=20.5 s/m, v c0 =6.2A m/s and T = 2.45 K we get E = 157 K. We 
thus have the important result that the energy barrier is reduced by the flow 
velocity in a parabolic form. We emphasize that this is the consequence of 
the linear variation of v c (T) and the temperature independence of the slope 
b of the slip rate. If b becomes temperature dependent the energy barrier as 
a function of velocity will contain higher order terms. 

We are now in the position to extend the analysis of the dc flow exper- 
iments of the Berkeley group (Fig. 10 of Ref. ^). In these experiments the 
energy barrier E* is obtained at various constant temperatures for a certain 
range of velocities, where E* is defined as 

E*(v,T) = E(v) — T ■ lnT. 



E(v) is the kinetic energy barrier and T is an "attempt frequency" (in Hz), 
see Fig. ||. Using the same energy scale as with the slip rate, namely E(v) — 
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Fig. 3. Results for the energy barrier E* versus dc flow velocity at vari- 
ous temperatures ranging from 1.83 K (lower left) to 0.35 K (upper right) 
obtained by the Berkeley group (digitized data from Fig. 10 of Ref. ||). 

E(v c (T)), we calculate the energy difference AE* between two isotherms T\ 
and T 2 at the same velocity v to be 

AE* = -£(« c (Ti)) + E(v c (T 2 )) - {T x - T 2 )lnT 

where V is assumed to be constant in temperature. Inserting v c (T) into 
Eq. © gives E(v c (T)) = E T 2 /T$ and hence 

AE* = -Eq{TI - if )/2# - (Tx - T 2 ) lnT 

Introducing AT = T 2 — T\ as a new variable leads to 

AE* = -^{ATf + {^ + lnT)AT. (4) 

We construct now the universal E(v) curve for all dc data by connecting 
succesively all measured isotherms in Fig. ||| by a vertical shift to the data at 
the reference isotherm T 2 = 1.83 K, see Fig. |j. The resulting curve follows 
approximately the quadratic variation of Eq. (^), but a small third order 
term (15% at the highest E* values) is needed for a perfect fit. In Fig.|| the 
vertical shifts AE* are shown and Eq. @ is fitted. This yields Eo = 57 K 
(for T = 2.83 K extrapolated from Fig. 9 of Ref.|) and T = 1.2 • 10 4 Hz. 
Although these values differ from those of the ac experiments discussed above 
the general properties of both experimental results are similar. 
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Fig. 4. Connecting all isotherms of Fig.|| by vertical shifts produces an 
universal velocity dependence of the energy barrier. The fitted curve is a 
third order polynomial, see text. Energy zero chosen at v/v c q = 1. 
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Fig. 5. Energy shifts AE* obtained from constructing the connected curve 
of Fig. | versus the temperature difference from 1.83 K. Equation (4) is fitted 
to the data. 
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4. Conclusion 

Starting from the correct relation between cumulative distribution func- 
tion and slip rate we have re-analyzed the published data on single phase 
slip events. We find a slip rate which varies exponentially with flow velocity. 
We deduce the energy barrier for vortex nucleation, decreasing quadratically 
with flow velocity. Finally, we extend the analysis of the dc flow experiments 
and obtain a consistent picture for the results of both methods. We do not 
offer any physical model for the results of our analysis. Further work on 
vortex nucleation is needed to explain the linear temperature dependence of 
the critical velocity, the exponential increase of the slip rate or the quadratic 
dependence of the energy barrier on the flow velocity. 

Returning to our own experiment on the transition to turbulence in 
the oscillatory flow of superfluid 4 He at millikelvin temperaturesS we note 
that in our case the failure rate (in units of s/m) of the laminar phase is 
proportional to the velocity increase above the critical velocity for turbulent 
flow. The lifetime (in units of s), however, may become metastable being 
limited ultimately only by creation of vorticity due to natural background 
radioactivity or cosmic rays. The lifetime of the turbulent phase grows 
rapidly with increasing drive power, becoming infinite at some critical value. 
The stability of both the laminar and the turbulent phases can probably not 
be understood in terms of single vortex nucleation effects but rather by the 
dynamics of fully developed superfluid turbulence. 
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